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We briefly review and compare three methods (one perturbative, one based on Ward Identities and one non- 
perturbative) for the calculation of the renormalization constants of lattice operators. The following results 
are presented: (a) non perturbative renormalization of the operators with light quarks; (b) the renormalization 
constants with a heavy (charm) quark mass and its KLM improvement; (c) the non perturbative determination 
of the mixing of the AS = 2 operator. 



1. Introduction 

The Renormalization Constants (RC's) of lat- 
tice operators are a necessary ingredient in the 
calculation of physical Matrix Elements (ME) 
from lattice correlation functions. Schematically, 
a physical amplitude A a ^p of a transition a — > (3 
is given by 

A a ^ p = C w {fi)Z {afi)(a\0{a)\(3) (1) 

Cw is the Wilson coefficient of the OPE, calcu- 
lated in Perturbation Theory (PT) at a renormal- 
ization scale /i; (a\0\(3) is the lattice bare ME of 
the operator relevant to the physical process, sim- 
ulated numerically; a is the lattice spacing; Zq 
is the operator's renormalization constant. It can 
be calculated either in PT, or, non perturbatively, 
through Ward Identities (WI) or with a recently 
proposed non-perturbative (NP) method (Sects. 
2-4). This review is focused on the systematic 
errors which afflict these calculations (Sects. 5 - 
8) It is limited to Wilson and Clover 2- and 4- 
fermion operators. 

2. Evaluation of RC's from PT 

The perturbative renormalization of lattice 
operators will be reviewed elsewhere [1]; here 
we only highlight the essential features of the 
method. Consider a generic bilinear bare 
fermionic operator 0(a) = ipTip and a generic 1PI 
Green function To(p, go, a) = (p|Or|p) (go is the 
bare coupling and p stands for the momenta of 
the external states). We will denote by O(ji) and 



fo(p, g 2 ,fi) the finite renormalized counterparts 
of the above quantities and by Tq^ the tree level 
estimate of To- The RC Zo(a,yL) is then defined 
through f (p, g 2 ,fi) = Z {gl,a^)T {p,gl,a). It 
depends on the choice of regularization and renor- 
malization scheme. The 1-loop estimate of To is 

r ° = [1 + (Jp (Tln(pa)2 + Clat Wo ] (2) 

where 7 is the anomalous dimension and Clat 
a finite coefficient. The renormalization condi- 
tion usually imposed on To m lattice perturba- 
tive calculations (e.g. [2]) is the one chosen in 
the corresponding perturbative calculation in the 
continuum; normally this is the MS renormal- 
ization scheme in a Dimensional Regularization 
(DR), such as Dimensional Reduction, 't Hooft - 
Veltman etc. Thus we impose 

2 

lim a ^ t o \ v 2 =li 2 = [1 + JL^Cdr]^ (3) 

where Cjjr is the finite coefficient of the 1-loop 
calculation in the continuum (i.e. analogous to 
Clat in eq.(2)). This is a choice of convenience: 
the Wilson coefficients Cw of eq.(l) are known in 
MS (e.g. [3]). The RC's must consequently be 
in MS if the l.h.s. of eq.(l) is to be Renormal- 
ization Group invariant. The 1-loop perturbative 
estimate of the RC then follows: 
2 

z = i- [7 MH 2 + C LAT - C DR ] (4) 

The dependence of Zq on Cjjr comes from 
the choice of the MS renormalization scheme, 



whereas its dependence on Clat from the choice 
of lattice regularization. Two perturbative calcu- 
lations are thus necessary, one in the continuum 
for Cjjr and one on the lattice for Clat- 

The following points must be stressed: (1) The 
reliability of the perturbative calculation rests on 
the assumption that we work with a large enough 
cutoff; i.e. Aqcd "C a -1 . (2) The renormal- 
ization condition of eq.(3) is defined in the limit 
o — > 0; i.e. all 0(a) terms (and higher) vanish. 
(3) Operators which are protected by WI's, such 
as the vector current V M = tpj^tp, have 7 = and 
Cjjr = (for the axial current = ipj^jsip, 
a finite renormalization is neccessary in order to 
impose the WI's). In this case the MS renormal- 
ization condition of eq.(3) coincides with the one 
used in the NP method (see Sect. 4). 

The calculation of the RC's in 1-loop PT suf- 
fers from 0(<7o) systematic errors. Moreover, the 
lattice coupling g$ is known to be a bad expan- 
sion parameter [4]. One particularly simple rem- 
edy consists in the choice of the boosted coupling 
= gg/Un (Un stands for the average plaque- 
tte) as a better expansion parameter. This we 
call Boosted Perturbation Theory (BPT). Note 
that in theory, when working at 1-loop, it is im- 
material which coupling (e.g. 9o,9jj-g, 9o) is used 
in eq.(4). This ambiguity in the choice of coupling 
is lifted in (and up to) 2-loop PT [5]. 

Another subtlety is that Standard lattice PT 
(SPT) ignores finite volume effects, since loop mo- 
menta are integrated rather than summed. Re- 
cently [6] this was taken into account by summing 
over discrete momenta; this is known as Discrete 
SPT (DSPT). Combining BPT with Discrete PT 
gives the RC's in Discrete Boosted PT (DBPT). 

The PT method is reliable for finite and loga- 
rithmically divergent operators, the renormaliza- 
tion of which does not require a subtraction of 
lower dimensional operators. In the case of mix- 
ing with operators with lower dimension, a non- 
perturbative subtraction is necessary. This is be- 
cause any non perturbative contribution to the 
mixing coefficients, being ~ exp[— J 90 dg/(3(g)], 
when multiplied by the divergent behaviour o _1 , 
will give rise to an o _1 exp[— J 90 dg/(3(g)] ~ 
Aqcd term, which will not vanish in the con- 



tinuum limit o — > 0. This is the most serious 
limitation of the perturbative method. 

3. Evaluation of RC's from chiral WI's 

The most accurate method for a non pertur- 
bative determination of RC's is based on chi- 
ral WI's. Only a limited number of RC's can 
be extracted from them, such as Zy , Za and 
the ratio Zp/Zs of the RC's of the pseudoscalar 
(P = ipj 5 ip) and scalar (5 = tptp) densities. The 
RC's are determined by requiring that bare cor- 
relation functions calculated (numerically) on the 
lattice, when renormalized, obey the WI's of the 
theory. For example, the WI related to the vec- 
tor flavour symmetry on the lattice gives rise to 
a point split conserved current V£ [7], for which 
Zyc = 1. Since the local vector current is not 
conserved on the lattice, we fix its renormaliza- 
tion by requiring that the renormalized current 
Vft = ZyVp obey the same vector WI. This im- 
plies that Zy can be determined from a ratio of 
ME's of the kind 



(5) 



For the axial current and the ratio P/S a sim- 
ilar procedure can be applied, provided that the 
breaking of axial symmetry by Wilson fermions 
is properly restored at vanishing lattice cutoff 
[7]. If the calculation is performed with WI's on 
hadronic states, everything is explicitly gauge in- 
variant. If the external states are quark states, 
gauge fixing is necessary and the presence of 
Gribov copies may in principle undermine the 
method. 

4. Non-perturbative evaluation of RC's 

Recently, the following new non-perturbative 
(NP) method, based on the renormalization of 
quark state correlation functions has been pro- 
posed [6]: the quark state correlation function 
G {x,y) = (0\xj}(x)O r {0)^(y)\0) is calculated nu- 
merically in coordinate space at fixed coupling 
(3 and hopping parameter K. Its forward ampu- 
tated Fourier Transform, Ao(p), is then calcu- 
lated. From it, by a suitable projection in Dirac 



Method 


Z V {W) - Ref. 


Z V {SW) - Ref. 


SPT 
BPT 
WI(2-pt) 
WI(3-pt) 


0.83 [2] 
0.71 [4] 
0.57(2) [8] 
0.72(3) [8] 


0.90 [11] 
0.83 [4] 
0.82(1) [9] 
0.82(1) [9] 



Table 1 

Z v calculated in SPT, BPT and with the WI 
method from ratios of (V M V M ) (2-pt) and (PV M P) 
(3-pt) correlation functions with the Clover(SW) 
and Wilson(W) actions. 



indices, P, we obtain the projected correlation 
function To = Tr[PAo], on which the following 
(off-shell) renormalization condition is imposed 



Zo {afi, gDZ. 1 {afi, gl)T {pa) | p 2 =M 2 



(6) 



If we know Z^ we can solve eq.(6) for Zq- The 
best determination of Z^ is from the same eq.(6) 
with O = Vf? , for which Z v c = 1 (equivalently 
we can use O = V M and Zy from of eq.(5)). 
Another estimate of Z^ can be obtained from 
a lattice version of the propagator's derivative 
"dS^/dp" (see [6]). Note that the result de- 
pends on the momentum of the external states 
and on the gauge. The Landau gauge is a good 
choice. 

In the case of V M and A^, the renormalization 
condition of eq.(6) defines renormalized currents 
which automatically obey the Current Algebra 
(see [6] for the proof and some practical restric- 
tions). It should also be noted that Cjjr = for 
the ME's of the currents and thus the two renor- 
malization conditions eq.(3) and eq.(6) coincide. 
For all other operators, the renormalization con- 
ditions being different, we still need to match the 
RC computed from eq.(6) to the continuum RC of 
eq.(3). This matching has to be computed in con- 
tinuum PT. Thus with the NP method we avoid 
lattice PT with its problem of a bad expansion 
parameter. 

5. Discretization Errors 

At finite lattice spacing, operators mix with 
higher dimensional (irrelevant) operators. This 
mixing spoils the multiplicativity of the renor- 
malization of the lattice operators, which is only 



recovered in the limit o — > 0. The contribution 
of higher dimensional operators to the numerical 
calculation of RC's (from WI's or NP-ly) depends 
on the ME's used. This dependence signals the 
presence of systematic errors which are effectively 
0{ap) in the scaling limit. With the Clover ac- 
tion, the RC's of improved operators (in the spirit 
of [10]) have effectively 0{g^ap) errors. In order 
to control these systematic errors we must ensure 
that fi <C 0{a~ 1 ). Since we also need a good 
control of non perturbative effects, working in a 
window 



Aqcd <M < 0{a~ 1 ) 



(7) 



is an essential ingredient for the reliability of 
any renormalization procedure performed at fi- 
nite cutoff (e.g. BPT). 

In conclusion, when the condition of eq.(7) is 
satisfied, the estimates of the RC's from 1-loop 
PT (good up to 0{g^)) and those obtained from 
different ME's either with the NP or with the WI 
method (good to 0{ap) in the Wilson case and 
0{g^ap) in the Clover one) should all agree. In 
present day calculations at (3 = 6.0, however, such 
discrepancies are 20% - 30% for Wilson fermions. 
They reduce to 5% - 10% for light Clover fermions 
(see Table 1). Notice that BPT also improves the 
agreement. This good control of the systematic 
errors with the Clover action at current (3 values 
can worsen, however, as the quark mass increases. 
This case will be considered in detail in Sect. 7 

6. Light Quarks 

Results from WI estimates for the Wilson and 
Clover action for quark bilinear operators have 
been around for some time (see Table 1). Last 
year the first NP results have been presented 
with the Clover action (/3 = 6.0, K = 0.1425, 
V = 16 3 x 32 and 36 confs.) [6], [12]. The esti- 
mates of Z v and Z s obtained from BSPT, BDPT, 
NP and WI (the latter only for Zy) were in good 
agreement in a window of /i 2 o 2 values. The sit- 
uation was less satisfactory for Za, where lattice 
artifacts were seen. For Zp a discrepancy be- 
tween the NP estimate and those from BSPT and 
BDPT signals large 2-loop corrections. 

In this conference the first results for the same 
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Figure 1. Wilson action estimate of Zy from NP 
method (points), from PT (continuous line) and 
from tadpole improved PT (dashed line). 



operators at the same V and (3 values (K = 
0.1515, 0.1530, 0.1550, and 125 confs.) were pre- 
sented for the Wilson case [13]. A good agreement 
between the SPT and NP estimates is seen in a 
large window in [i 2 a 2 (Figs 1,2 at K = 0.1530). 
The Za results (not shown here) are of similar 
good quality, whereas Zs (also not shown) fluc- 
tuates more (just like in the Clover data of [6]). 
These Wilson results (but obtained on a subset 
of 36 confs.), compared to the Clover ones of [6], 
are significantly less fluctuating. This may be at- 
tributed to the fact that the Clover action and the 
improved operators are less local than the Wilson 
ones. 

The same collaboration has also calculated nu- 
cleoli structure functions [13], [14]. This requires 
the renormalization of operators like and 
0{iiva}-> from which the 1st and 2nd moments of 
the structure functions are obtained. For exam- 
ple, is defined by 



= V'T{m- D ^}V' - traces 



(8) 



where {. . .} stands for symmetrization of the 
Lorentz indices. The PT estimates of these RC's 
with the Wilson and Clover actions have been 



Figure 2. Wilson action NP estimate of Zp 



obtained in [15] and extended to higher rank op- 
erators in [14] for the Wilson case. 

The RC of {flu} (K = 0.1530) is given in Fig. 
3; similar results have been obtained for the other 
relevant operators. Again a satisfactory window 
in [i 2 a 2 is visible, in which the operator's anoma- 
lous dimension is clearly seen in the slope of the 
curve. 

So far we have assumed that the Gribov ambi- 
guity does not introduce any errors. This assump- 
tion may be questioned in the light of recent re- 
sults which demonstrated that: (1) at least some 
Gribov copies are lattice artifacts [16], [17] and (2) 
they spoil the measurement of some gauge de- 
pendent correlations [16]. In [18] the influence of 
Gribov copies was explicitly studied by evaluating 
Za, on the same ensemble of 36 confs. (/3 = 6.0, 
V = 16 3 x 32, K = 0.1425), both from a gauge 
independent WI on hadron states and a gauge de- 
pendent one on quark states. In the latter case 
Za was computed on several Gribov copies. The 
result 



Za = 1.06(2) gauge ind. 
Za = 1.08(5) gauge dep. 



(9) 



supports our assumption. We will assume that 
this result generalizes to all RC's calculated either 
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Figure 3. Same as Fig. 1 for the RC of the oper- 



ator O 



{/ii/} 



from quark state WI's or from the NP method. 

7. Heavy quarks 

As the quark mass increases towards the charm 
(e.g. in semileptonic decays) the systematic error 
of 0(am) in Wilson fermions (O^g^am) in Clover 
fermions) is bound to grow beyond control. In [4], 
[19] a remedy has been proposed in the form of 
the so-called KLM normalization factors. These 
arise when a tree level lattice object (propaga- 
tor, correlation function etc) is matched to its 
continuum counterpart at zero spatial momentum 
p = 0. For example consider the continuum tree 
level 2-quark correlation function of the local vec- 
tor current in 3-momentum space (p = 0): 



C. 



CON 



I*' I* 



vWiWMoWjMIo) (io) 



where 1 and 2 label distinct flavours. In terms of 
its lattice counterpart C LAT (p = 0,t), calculated 
with Wilson fermions, it is given by 

C% N = 2K X 2K 2 {1 + omi)(l + am 2 )C V LAT (11) 

Besides the traditional 2K quark field normaliza- 
tions, there are the extra (1 + am) KLM normal- 
izations. They differ from 1 at a ^ 0. The above 



relationship is also valid when we match C^ k ON 
with C^ AT , the quark correlation function of the 
spatial components of the conserved current V£ 
(k = 1,2,3). For the temporal case V C there is 
an extra factor, as observed by Bernard [20]: 



C. 



CON 



2Ki 2K 2 (1 + omi)(l + am 2 ) 



ami + am 2 



LAT 



(12) 



The claim of [4], [19] is that taking these KLM 
corrections into account when calculating, say, 
Zy from eq.(5), will soften its "mass dependence" 
(which is really a systematic 0(am) error). A 
counter - claim raised in [21] is that since KLM 
factors only remove the 0(am) error at tree level, 
they might at best cause a softer "mass depen- 
dence" of the RC's obtained at a given, fixed ME. 
They will not, however, eliminate the 0(am) sys- 
tematic error which shows up as dependence of 
the Zo estimate on the ME's it was obtained 
from. Also, since the KLM factors are calculated 
at p = 0, they do not correct C(o|p|) effects. 

In the results shown below [22], the value of the 
quark mass used in the KLM factors is simply 
taken to be am = 1/(2K) - 1/(2K C ) with K c 
extracted non perturbatively. Using, instead, the 
original Mean Field Tadpole Improved (MFTI) 
quark mass of [4] gives similar results [21], [23]. 

In Fig 4 we reproduce from [21] (/3 = 6.4, V = 
24 3 x 60, 20 confs.) the "mass dependence" of the 
ratio 



Ri 



(K\Vq C \D) 
(K\Vo L \D) 



(13) 



(obtained from 3-point correlation functions) 
without and with KLM corrections. It is clear 
that the KLM claim is well supported. Next we 
consider the ratios of 2- and 3-point correlations 



(Q\v k c \p) 

m L \ P ) 



and Rs 



(K\V k C \D) 
(K\V k L \D) 



(14) 



for which the KLM factors cancel out from nu- 
merator and denominator. Comparing the re- 
sults of R 2 and Rs (Fig. 5) to those of Ri (Fig. 
4), we see how the KLM prescription fails: The 
KLM corrected #i ~ 0.75 and R 2 ~ 0.65, al- 
though essentially independent of the quark mass, 
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Figure 4. Zy calculated from Ri without and 
with KLM corrections for Wilson fermions. 



Figure 5. Zy calculated from R 2 and R% for Wil- 
son fermions. 



give different estimates for Zy . The situation is 
even worse for the Zy estimate obtained from R3. 
Thus, even at (3 = 6.4, there are big 0(am) sys- 
tematic errors unaccounted for by tree level KLM 
factors. 

The same KLM ideas have been applied and 
tested in the Clover case. Here there are extra 
KLM factors arising from the Clover field rota- 
tions [10] V —> i> R = [l - fr 5]V>. As an 
example we give the KLM factors of the quark 
correlation function of the local improved current 
(01^00^(0)^)10); 
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(15) 



Different KLM factors arise in the case of the con- 
served improved current V^ 1 (see [22] for details). 

Since the KLM factors respect Clover improve- 
ment, they are 1 + 0(a 2 m 2 ). Thus, they can- 
not rectify the O(gQam) corrections present in the 
Clover approach. This is demonstrated in Figs. 
6, 7 (f3 = 6.0, V = 16 3 x 48, 30 confs.) where we 



plot the ratios 

(Q|^ ci Ip) 
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and R5 



(K\V k CI 
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(K\Vl\D) 



(16) 



without and with KLM corrections [22]. The 
large linear O^g^am) term of the KLM corrected 
Ri (Fig. 6) cannot be offset by the 1 + 0{a 2 m 2 ) 
KLM factor. In Fig. 7, the data does not fall 
on a smooth curve because R$ also depends on 
a 2 (mi — m-i) 2 (it is symmetric under mi <-> m^). 
Comparing the KLM corrected data of Figs 6 and 

7, we see that the common linear behaviour up 
to m ~ 0.5 is followed by a non universal, ME- 
dependent behaviour at large quark mass. This 
systematic error cannot be corrected by KLM fac- 
tors 

8. AS = 2 

The reliability of the evaluation of the Bk pa- 
rameter with Wilson fermions has been plagued 
by the long standing problem of the bad chiral 
behaviour of the ME (K \O A5= 2 |X°) . The non 
vanishing of the above ME in the chiral limit 
stems from two systematic effects. The first is 
the 0(<7o) error in the 1-loop perturbative cal- 
culation of the logarithmically diverging RC, Zq, 
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Figure 6. calculated from i?4 without and 

with KLM corrections for Clover fermions. 



Figure 7. calculated from R$ without and 

with KLM corrections for Clover fermions. 



and the finite constants, Zy. (k = 1,2,3), which 
determine the mixing of the Oas=2 operator with 
dimension-6 operators of the wrong naive chiral- 
ity: 

Oa5= 2 (m) = Z (a)[O AS= 2( a ) + ^i(a)0 5 p(a) 

+Z 2 (a)0 VA (a) + Z 3 (a)Os PT (a)} (17) 

(the definitions of the above operators can be 
found in [26]). The second systematic error comes 
from 0(am) discretization effects in the four ma- 
trix elements (K \O AS= 2 \K°) and (K°\O k \K°) 
(k = SP, VA, SPT), which are calculated nu- 
merically at finite lattice spacing. In a pre- 
liminary comparative study of Bk with Wilson 
and Clover fermions [24], [25], the latter source 
of error was investigated. The implementation 
of the Clover action caused a very marginal im- 
provement of the chiral behaviour of the ME. 
Thus the principal source of error seems to be 
the perturbative determination of the RC's. In 
[26], [27] the NP method was applied to this op- 
erator (for details see [26]): The amputated 4- 
point quark correlation functions of the opera- 
tors of eq.(17), suitably projected with projectors 
Pas=2 and Pj(j = 1,2,3) are calculated numer- 
ically. By imposing the following four renormal- 
ization conditions on the subtracted correlation 
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Tr[PjA s ] = 0;j 



ZAS=2Zj 2 Tr[P AS= 2A s 



(18) 



we obtain the four Z's. In Fig. 8 we report from 
[26], [27] a preliminary but very encouraging re- 
sult (at (3 = 6.0 and V = 16 3 x 32) with Clover 
fermions and NP mixing constants. It is clear 
that the NP mixing reproduces the expected chi- 
ral behaviour of the ME, whereas SPT and BPT 
do not. This is a crucial test that Wilson fermions 
have to pass if they are to predict the Bk parame- 
ter with the same accuracy as staggered fermions. 

9. Conclusions 

For light quarks, a number of improvements 
(BPT, the Clover action, the NP method) have 
brought the different estimates of the RC's to 
close agreement. The newest encouraging results 
in this area come from the NP calculation of the 
RC's of the operators which determine the nu- 
cleoli structure functions. 

The NP estimates of the mixing coefficients of 
the Oa5=2 operator apparently reproduce the ex- 
pected chiral behaviour of the ME. If this pre- 
liminary result is verified by the more systematic 
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Figure 8. Chiral behaviour of (0+) = 
(K \O AS= 2\K°)/(P) 2 with mixing constants de- 
termined from SPT, BPT and NP methods. f K 
and Mk are the decay constant and mass of K°. 



runs currently performed by the APE group, it 
will reopen the way to the long awaited A J =1/2 
calculation. 

For charmed quark masses, WI estimates of 
RC's suffer from systematic discretization errors 
at the current (3 values. These are only partially 
accounted for by tree level matching. 
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